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Abstract. - 

New developments have spurred interest in magnetic moments (/i-s) of baryons. The mea- 
surement of some of the decuplet fi-s and the findings of new sumrules from various methods 
are partly responsible for this renewed interest. Our model, inspired by large colour approxima- 
tion, is a relativistic self consistent mean field description with a modified Richardson potential 
and is used to describe the fi-s and masses of all baryons with up (u), down (d) and strange 
(s) quarks. We have also checked the validity of the Franklin sumrule (referred to as CGSR 
in the literature) and sumrules of Luty, March-Russell and White. We found that our result 
for sumrules matches better with experiment than the non-relativistic quark model prediction. 
We have also seen that quark magnetic moments depend on the baryon in which they belong 
while the naive quark model expects them to be constant. 
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Introduction. 't Hooft suggested that the inverse of number of colours (N c ) could 
be used as an expansion coefficient in the otherwise parameter free QCD [1]. Based on 
this, Witten [2] suggested mean field (MF) description for baryons - prompting the use of 
phenomenological interquark potential tested in meson sector. Indeed baryon mass (M) was 
calculated at MF level [3] using Richardson potential as an interquark one [4] . The potential 
has confinement and asymptotic freedom (AF) built in - with a single scale parameter (A ~ 
400 MeV). Although the lattice QCD gives a confinement scale (~ 400 MeV) - the AF 
scale as given by perturbative QCD is ~ 100 MeV . For computing the hadron properties 
aA~ 400 MeV is required whereas for high density strange quark matter (SQM in short) 
A ^ 100 MeV [5]. This is not surprising as for the SQM the confinement gets screened and a 
much reduced A appropritate for AF part is required. Hence it is important to separate out 
the two scales and re-do the hadron properties -like magnetic moments and masses for the 
baryons. In fact, Bagchi et al [6] calculated the masses (M) and magnetic moments ^ of A ++ 
and f2 _ with a modified two parameter Richardson potential using relativistic Hartree-Fock 
(RHF) method. 

The known /i-s cannot be explained by any model exactly and in addition there are some 
sumrules which are interesting to look at. The first one that is much referred is the Franklin 
sumrule, named erroneously after Coleman and Glashow [7]. The latter authors derived mass 
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sumrules for SU(3) breaking and found specific /U-s but not the actual sumrule for the octet 
/i-s - which was first done by Franklin [8] and so should carry his name( 1 ). 

Presently, extending the calculation of Bagchi et al [6] with the modified potential we find 
all baryonic (decuplet and octet) /i-s and check few sumrules to test symmetry assumptions. 
We compare with experiment and other calculations hoping more experimental decuplet /i-s 
to be soon deduced. Our values are thus predictions. One can view our calculation in the 
spirit of large N c or consider it as a relativistic MF calculation with a potential having AF 
and confinement property built into it. We compare our results with those obtained from 
analytic large N c formulations which are exact calculations but assumes static limit as M is 
infinite for N c = oo. Our Ms are finite, so the comparison is interesting - displaying the 
effect of the wave functions and whatever dynamics it contains. 

Formalism. - One needs to sum all planar gluon exchange diagrams to deduce an ef- 
fective interquark potential. Analytic derivation of such a sum being absent, a potential like 
Richardson potential is chosen from meson phenomenology and then tested for baryons and 
quark stars [3, 5] 

V (r) = - ^yjl- 67 L r [A 2 r-/(Ar)l (1) 
v ; 2N C 33 - 2N f w 



jVc + l 

2N C 



is colour contribution, N c is number of colours (= 3), Nf is number of flavours (= 3). 

f( t) -i_4 r dq ex p(-^ (2 ) 

Richardson calculated, non-relativistically, the masses of two heavy mesons, J /iff and T [4]. 
Crater and Van Alstine [9] obtained masses of both light and heavy mesons using a relativistic 
two body Dirac equation with A = 401 MeV. Dey et. al. [3] calculated the baryonic 
properties like M and \i of fl~ with A = 400 MeV using RHF method. However, when the 
potential was used in strange star calculation [5] the required A ~ 100 MeV. Strange stars 
are very compact objects composed of high density SQM. Debye screening length of the gluon 
suppresses the confinement due to the medium effect and a lower A is sufficient. 

Separating the confinement and AF scales the modified potential is given by 

A' = 350 MeV (confinement) and A = 100 MeV (AF) gave satisfactory results for the 
hadrons A ++ and fl~ [6]. Also the strange star calculation improves [10], - increasing the 
value of the strong coupling constant a s in Debye screening term. This is more consistent 
in the framework of the star calculation in view of the findings of [11] where the density 
dependence of quark masses is found from a s . 

We set out to apply this potential to find the /z-s of all other baryons, and check magnetic 
moment sumrules, with no more extra parameter to adjust. 

Details of Calculation. - The Quark wave function for quarks in lowest (lsi/2) orbital 
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u,d,s (4) 



( 1 )We are grateful to Jerry Franklin for pointing this out to us and allowing us to correct the mistake. 
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Xm is Pauli spinor, a is Pauli matrix. The Hamiltonian and corresponding HF equations are: 

3 

i i<j 

z = 1,2,3; j = 1,2,3; U = a.p i J rf3m l (6) 

[t q + UJ q (n )}(f) q (ri) = £ q (j> q {r X ) (7) 

ti is the kinetic energy operator for the ith quark, V(j*ij) is the modified potential, e q -s are 
the single particle energies and uj q -s are the single particle potentials. 

Wqi(ri) = / 4> q 2*{r2)V{r 12 )4> q 2(r2)r2 2 dr2 + <j> q 3H r 3)V(r 13 )<f) q3 (r 3 )r 3 2 dr 3 (8) 



Using the wave functions given in eqn. 0] , we get sets of coupled differential equations 

dG 

—A - (m q -ut q + eg) F q = 
d - + ( - ; F <? + ( £ ? - ^ - m ?) G ? = (9) 



From e q -s, the energy is obtained from the following equation 

E = e u + e d + e s - i J (j)^ (ri)uj u (j) u {ri)rldri - 

^ y (r 2 )uj d <j) d (r2)rldr2 - i J fa* (r 3 )u a <l> a (r 3 )r§dr 3 (10) 

Here the subtractions are to avoid double counting which comes when applying variational 
principle to derive HF equations. 

The coupled differential equations are to be solved self consistently. Procedure is to check 
the convergence in the energy value. Moreover, to make the calculation more transparent 
and easier we take recourse to expansion of wave functions in oscillators and subsequent 
diagonalisation, 

G(r) = J2°n R n0 (H) 

n 

F{r) = J2 D ^Rrni (12) 

m 

C-s and D-s are coefficients. This reduces the differential equations to an eigenvalue problem. 
Starting with a trial set of C-s and D-s, the solution is found self consistently by diagonalizing 
the matrix and putting back the coefficients till convergence is reached. In general R n i (r) is 



/?,„(/) = J T{ J + n i +3 _/ exp(-^)L l n + Hr^ (13) 



L n 2 (r 2 ) are Associated Laguerre polynomials. In the calculation, r is replaced by r/b, b 
being the oscillator length which may be different for G(r) and F(r) - b and 6', respectively. 
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The centre-of-mass (CM) momentum is not well-defined in RHF solutions and this entails 
a spurious contribution from the CM kinetic energy to the total energy. Since the relative 
importance of this effect increases as the number of particles decreases, it is necessary to 
correct it for systems formed of few particles. This has been done here by extending the 
Peierls-Yoccoz procedure of nuclear physics. The spurious contribution is denoted by Tcm 
and the baryon mass M is E - Tcm- The Peierls-Yoccoz procedure corrects the energy but 
for /i - correct boosted wave functions are needed. This is discussed in the reference [12]. 
Boosting is not attempted in the present paper since our ultimate aim is to use the procedure 
for a large system, a massive star, where these corrections are not relevant. Therefore, our 
results on masses of decuplet and octet baryons, take care of center of mass correction but 
not for the magnetic moments. Still, we are able to compare our calculated fx — s with infinite 
mass static large N c models which also do not involve centre of mass correction( 2 ). 

Energy differences between decuplet and octet baryons are obtained using the simplest 
idea of instanton physics [14] where a is the instanton induced potential between a (u,d) pair 
and [3 is the same between a (u, s) or (d, s) pair. 



E N = E a - 3a/2 (14) 

E A = E A -a-P/2 (15) 

E s = Ee. ~ 3/3/2 (16) 

E E = Ew* - 3/3/2 (17) 



There are alternative methods also e.g. incorporating colour magnetic interaction energy [15]. 
The r.m.s. radii r av can also be estimated using the expression 




" [(G ll (r) 2 + F ll (r) 2 ) + (G d (r) 2 + J F rf (r) 2 ) + (G s (r) 2 +F s (r) 2 )]r 4 dr (18) 



But as the potential is spin independent, the radii for octet and decuplet members become 
the same. 

In quark model, the magnetic moment associated with a quark is given by [16]: 

^ = y / ( Fx j) d * r (is) 

e q is the charge of the quark, j is the current associated with the quark. fj, q reduces to : 

2 r°° 

Mf (1)] = - (+K 3 J G(r)F( r ydr (20) 

Baryonic /i is found using baryonic wave functions ^ S pin X ^flavor [17]. 

Results. Subtracting Tcm from HF energy E, we obtained the baryonic masses M 
which are given in table [i] along with the experimental values. Tcm lies ~ 100 MeV for all 
the baryons. We have chosen a — 188 and (3 = 103 MeV [14] which gives the overall 
best fit. The values of the oscillator parameters b and b' are chosen such that E becomes 
independent of variation of 6, b' . There are no other free parameters to fit. We have also 
found that r av is ~ fm and it decreases with increasing M . 



( 2 )Sincc an oscillator basis is used for the calculation, it is possible to do standard CM correction exemplified 
by Elliott and Skyrme [13] . 
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Table I - Charge averaged baryon masses (both decuplet and octet members) where A' is 350 MeV and 
A is 100 MeV; using 7x7 matrices. The quark masses are m u ,m,d ~ 10 MeV , m 3 = 150 MeV. 
We have chosen a — 188 MeV and (3 — 103 MeV for a good fit of nucleon and E comparing to the 
experimental value. 



Baryons 


Experimental Mass 
(MeV) 


Theoretical Mass 
(MeV) 


b 


b' 


A's 


1232 


1251 


0.83 


0.60 




1383 








E*'s 


1384 


1361 


0.83 


0.60 




1387 








E*'s 


1532 


1455 


0.77 


0.60 




1535 








Q~ 


1672 


1556 


0.70 


0.60 






938 






N 


939 


938 


0.83 


0.60 




1189 








E's 


1193 


1188 


0.83 


0.60 




1197 








A 


1116 


1098 


0.83 


0.60 




1315 








H's 


1321 


1282 


0.77 


0.60 



Table [H] shows a comparison between our /i-s with experimental and other theoretical 
values. The agreement of Ms and /i-s from our result with those from experiments is not too 
unreasonable. Results obtained in QCD sum rule (QCDSR) approach are taken from [18] for 
the decuplet and from [19] for the octet. Results obtained by another large N c approximation 
are taken from [20] where the authors fit the octet and the £l~ /t-s to predict the other 
decuplet /x-s. Experimental values are taken from [21]. We have also compared our result 
with lattice [22] and chiral perturbation theory, \pt [23]. Dai et. al [24] fitted the octet /i-s, 
in their fit A and fit B, by adjusting 10 parameters and then predicted the unknown /i-s. For 
the decuplet the agreement between the calculations and the three known /i-s are the only 
guides. 

The new experimental value of /ia++ agrees better with our result, QCDSR and lattice 
than the others. The Franklin sumrule [8] mentioned in the introduction is as follows : 

(/ip - Mn) + (ME- - + (MHO - MH-) = ^Franklin = 0. (21) 

The value of ^-Franklin is +0.48 from experimental /i-s, +0.47 from CDM [25] /i-s, +0.14 and 
+0.30 from Franklin's recent calculation [26] and +0.45 from our results. 

Taking infinite colour limit and the consequent static infinite mass one can get some other 
sumrules among the /i-s [27,28]. Luty et al. [28] coupled the mass expansion for the s quark to 
this limit which makes it interesting to compare with our calculation (tabic ITTTf t . Our results 
seem to be intermediate to the sumrules and experiment. 

With help of equation (|20|l we found that /i for a particular quark changes in different 
baryons, whereas quark /i-s are constant in naive quark model; fi u = + 1.852, fi u = — 0.972 
and fi s = —0.613 [21]. We found that with increasing M, the magnitude of quark /i decreases 
(see table Hv|. 
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Table II - Comparison of baryon magnetic moments found in different approaches. 





A++ 


A+ 


A u 


A- 


£*+ 


E* u 






H* 


Qr 


Ours 


5.77 


2.88 


0. 


-2.86 


2.81 


0.17 


-2.46 


+0.30 


-2.17 


-1.92 


Expt. [21] 


6.14 


2.70 


- 


- 


- 


- 


- 


- 


- 


-2.02 


QCDSR [18] 


6.14 


3.02 


0.0 


-3.07 


1.90 


-0.07 


-2.03 


0.80 


-2.71 


-2.02 


lattice [22] 


6.09 


3.05 


0.0 


-3.05 


3.16 


0.33 


-2.5 


0.58 


-2.08 


-1.73 


XPt [23] 


4.0 


2.1 


-0.17 


-2.25 


2.0 


-0.07 


-2.2 


0.10 


-2.0 


- 


1/JVc [20] 




3.04 


0.0 


-3.04 


3.35 


+0.32 


-2.79 


0.64 


-2.36 




Dai fit A [24] 


5.84 


















-2.08 


Dai fit B [24] 


5.86 


















-2.06 




P 


n 




£+ 


S° , A 


E~ 








Ours 




2.88 


-1.91 




2.59 


0.83, -0.71 


-0.92 


-1.45 


-0.62 




Expt. [21] 




2.79 


-1.91 




2.46 


-0.61 


-1.16 


-1.25 


-0.65 




QCDSR [19] 




3.04 


-1.79 




2.73 


-0.50 


-1.26 


-1.32 


-0.93 




CDM [25] 




2.79 


-2.07 




2.47 


-0.71 


-1.01 


-1.52 


-0.61 




Dai fit A [24] 




2.84 


-1.87 




2.46 




-1.06 


-1.28 


-0.61 




Dai fit B [24] 




2.80 


-1.92 




2.46 




-1.23 


-1.26 


-0.63 





Conclusions. The property of all baryons has been investigated with an improved 
Richardson potential in a tree level calculation in the large N c spirit and the results agree 
reasonably with experiments and other theoretical ones. The calculation is simple but it gives 
some of the quark dynamics which is absent in the exact calculation in infinite N c (infinite 
mass) static baryon model. Accurate determination of other decuplet fi-s will be helpful for 
testing models including ours which can be used to predict the equation of state for SQM. 

A different approximate relativistic many body method is perhaps possible for three quark 
systems but would not be relevant for us. We aim to test the validity of MF approach suggested 
by Witten [2] "QCD simplifies as N becomes large, and there exists a systematic expansion 
in powers of 1/N. In various ways, to be discussed later, this expansion is reminiscent of 
known phenomenology of hadron physics, indicating that an expansion in powers of 1/N may 
be a good approximation at 1/N = 3." and "The large N limit is, instead, given by a sort 
of Hartree approximation. The logic behind this approximation is as follows. For large N 
the interaction between any given pair of quarks is negligible - of order 1/N. But the total 
potential experienced by any one quark is of order one, since any quark interacts with N other 
quarks, each with strength 1/N. Thus, the total potential experienced by any one quark is 

Table III - Checking of six large N c analytic sumrules given by eqn (26) and (30) in Luty, March- 
Russell and White [28] with our results. 



Large N c analytic relations 


Our results 


Experimental 


(i) flp + \l n + Ms- = 


+0.05 


-0.28 


(ii) n s o — 2/i s - = 


-0.21 


+0.05 


(iii) fj, s - + fi s + - - ii v = 


+0.01 


+0.18 


( iv ) Hn- - Ms" - Ms- = 


+0.15 


-0.12 


(v) fl B a + 2/i E+ + 2/i E - + fj, n = 


-0.02 


-0.56 


(vi) jj, n - + 4^ s o — S)i s - + 8/i E + + 5fi s - - 3fi p + fi„ = 


-0.29 


-1.486 
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Table IV - Quark magnetic moments in our calculation 





A++ 


A+ 


A u 


A" 




E — 




E*~ Q~ 


tin 


+1.92 


+1.92 


+1.91 




+1.76 




1.64 




Md 




-0.96 


-0.95 


-0.95 




-0.87 




-0.82 


Us 










-0.72 


-0.72 


-0.67 


-0.67 -0.64 



of order one, but is a sum of many small, separately insignificant terms. As in statistical 
mechanics, when a quantity is a sum of many insignificant terms, the fluctuation around 
the mean value are very small. Thus, the potential experienced by one quark, apart from 
being of order one, can be regarded as a background, c-numbcr potential-the fluctuations are 
negligible. To find the ground state baryon, each quark should be placed in the ground state of 
the average potential that it experiences. By symmetry, the average potential is the same for 
each quark, so we should place each quark in the same ground state of the average potential". 
The beauty of MF approximation is that it can be used in both 3 quark system (baryon) and 
many quark system (quark stars). This has inspired our work. 
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